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Abstract 

Parity violation in a chiral, four-atom molecule is discussed. Given the geometrical positions 
of the four atoms, we calculate the anapole moment of it. This problem was first discussed 
by Khriplovich and Pospelov [lj. We give a detailed derivation for it so that it can be more 
accessible to wider range of scientists. We correct errors in their results and generalize their 
initial state to |s!/ 2 ) and \pi/f) states. We also discuss realistic candidates of the chiral molecules 
to which this approach can be applied. 

PACS numbers: 31.30.-i, 31.30.J- 

1 Introduction 

Charge-conjugation (C), parity (P), and time-reversal (T) violations in atoms and molecules are 
very attractive targets in searching for new physics beyond the standard model (BSM physics) as 
well as for their own developments of atomic and molecular physics and chemistry. For molecules, 
we can utilize selectively the precise spectroscopies and special environments of molecules suitable 
to determine fundamental properties of particles. A typical example is the recent improvement of 
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the upper bound on the electron electric dipole moment (EDM) by using a ThO molecule [2] ■ Also 
the measurements of energy difference between left- and right-handed molecules are on-going [3]. 
Such situations conventionally enforce on us complicated numerical calculations based on ab-initio 
molecular orbital (MO) calculation methods. Since there are so many papers on ab-initio MO 
methods, we cite only a textbook [3] and a recent paper [5J. Though they are crucially important, 
it is rather difficult to elucidate the physical origin of symmetry violations in a given molecular 
system from molecular orbital calculations. In order to remedy this deficit, analytical treatments 
of the same objects as those studied by numerical calculations have been awaited even with a cost 
of less rigor. From this point of view, the work [Ij by Khriplovich and Pospelov is methodologically 
very interesting. Given the geometrical positions of four atoms, they calculated the anapole moment 
of the chiral molecule which consists of the four atoms. There the P-violating (PV) anapole moment 
is directly and quantitatively related with geometrical structure of molecules. 

In this note we give a more detailed derivation of it than in the original work p], so that it 
can be more accessible to many physicists and chemists. Their final result is also corrected. 

We note here that in general, there are two possibilities concerning the origin of parity violation 
in molecules. One of the possibilities is the weak interactions which intrinsically violate parity. This 
effect is known to work in nuclei in atomic systems, which is discussed, for example, in Refs. mm- 
The other possibility is parity violation due to the geometric configuration of atoms in a molecule. 
When a molecule consists of more than four atoms, the molecule is not necessarily superimposable 
on its mirror image. Such a non-trivial transformation property under parity is called the chirality, 
and the chiral molecule and its mirror image are called enantiomers. When we take one of the 
enantiomers, parity is not a good symmetry to describe the quantum state of the chiral molecule 
since parity is violated by the configuration of the atoms. It is this latter case of parity violation 
which is discussed in Ref. [1! and which we consider in this paper. 

This paper is organized as follows. In Section 2, we describe our set-up of a four-atom molecule. 
In Section 3, we discuss the perturbation due to the Coulomb interactions between the valence 
electron and the light two atoms in the molecule. In Section 4, we calculate the anapole moment 
by using the results obtained in Section 3. We summarize our results and give discussions in Section 
5. In Appendix we discuss realistic candidates of the chiral molecules to which this approach can 
be applied. We work in the notation of the Landau-Lifshitz textbooks mm unless noted otherwise. 
We use the natural units in which h = c = 1 throughout this paper. 

2 Model Set-up 

Chirality in molecules first appears in four-atom molecules. Khriplovich and Pospelov considered 
a chiral molecule composed of four atoms, A\, A 2 , A 3 , A 4 , whose geometrical configuration is 
given in Fig. |T] without specifying its origin [lj. An assumption in this set-up is that Ai and A 2 
(whose electric charges are Z\ and Z 2 in units of the positron charge e, respectively) are light in 
comparison with A 4 and A 3 . Throughout this paper, we take the position of A 4 as the origin of 
our coordinate system, and the direction of F 3 as the 2 -axis. 

We first consider a valence electron which is captured by the diatomic molecule A 3 A 4 . To 
specify the position of the valence electron, we define the vector r whose initial and final points are 
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A 4 and the electron, respectively. We assume that the electron is first in one of the p 3/2 states with 
/i = ±3/2, where p is the z component of the total angular momentum. We also assume that the 
degeneracy in the energy levels of the atom A 4 is completely resolved by the presence of the atom 
A 3 , as assumed in Ref. |Tj. We treat the effect from the atoms A± and A 2 on the valence electron 
as a perturbation due to the Coulomb potential, 


V(r) = 


Z\ol Z 2 OL 
r — fi\ \r — 


(i) 


Then, the Coulomb interactions between the valence electron depicted as e in Fig. |T] and the atoms 
A\ and A 2 can be treated as a perturbation to the electron terms of the diatomic molecule A 3 A 4 , 
which induces a PV anapole moment. The unpaired electron is localized in the vicinity of A 4 , and 
hence its orbital can be well described in terms of the atomic orbitals of the atom A 4 . To lowest 
order of the approximation, the angular momentum about the A 3 A 4 axis is conserved. 



Figure 1: The configuration of the atoms in the four-atom molecule we consider in this paper. The 
numbers j (j = 1,..., 4) in the figure are the labels for the corresponding atoms Aj. e is a valence 
electron. The two planes are spanned by the atoms A 1 A 3 A 4 and the atoms A 2 A 3 A 4 , respectively. 
Also shown is the dihedral angle x °f the molecule. 

The final result we obtain in this paper is 


(a) = —jz [ri(s,p) - 4 r(s,p)\ ^2 Z X Z 2 
m e 45 


E s E p 


x (J • n 3 )(ni • [n 2 x n 3 ])[Ci(ri)C 2 {r 2 )fi 2 - C'i(r 2 )C' 2 (ri)ni] , (2) 

where a is the anapole moment operator and (a) its expectation value. Ry = o?m e /2 is the 
Rydberg energy (~ 13.6eV) with the electron mass m e and J is the total angular momentum, and 
the definitions of ni i 2 , 3 , Ci j2 , ri(s,p), r(s,p), E StP will be given in the next sections. The result is 
different from that of Ref. [lj in the factor 4 in front of r(s,p) and also in the overall factor by a 
factor of 2ao, where ao is the Bohr radius, ao = 1 /(m e a). Thus Eq. Q explicitly shows the relation 
between a PV observable and the geometrical structure of the molecule. In the subsequent sections 
we derive Eq. fl5J). 
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3 Perturbation due to Coulomb Interactions 


By the perturbation due to V(r), the states |p 3 / 2 ,^ = ±3/2) are slightly mixed with the states 
|sx/ 2 , d= 1/2) and |p!/ 2 , ±1/2). First we write the potential V(r) as 


V{r) = -Z x a\ 6 (r - n) ( r lm ± ^r lm n n -^—-^-- + 0(rl/r 3 ) 
1 ' dr m r 2 dr m dr n r 


n , \ i dll d d 1 3 3 

+ 0 (n - r) r m — - 1 - -r m r n — ---b 0 {r 3 /r\) 

dri m r i 2 dn m dr in n 


- Oi,Zi) -b (r 2 ,Z 2 ) , 


(3) 


where r* = |r)| (z = 1, 2) and the indices n and m (n,m = 1,..., 3) run over the three components 
of the Cartesian coordinates, ( x,y,z ). For later convenience, we integrate out the dependence of 
the potential on r (= |r|), and rewrite the potential Eq. (|3|) as 


V{n) 


Ry 



3 QmnlT'm'ri 



(4) 


where ft = r/r. The quantities D m and Qmn are defined as 

Dm — Z\C\{T\)Tl\m T j 

Qmn — ZiC2{ri)n lm n ln ~\~ ^2C2(j'2)'^ j 2m'^ j 2n ? 


coo 

C k (ri) = a 0 dr r 2 R 1 / 2 (r)R 3 , 2 {r) 
Jo 


y,h J,k 

' - d (r-r i ) + n —-9{r i -r) 


r»k- 1-1 


r k + l 


(5) 

( 6 ) 

(7) 


where Hi = r) /n (i = 1,2). The factors i? 3 / 2 (r) and i?i/ 2 (r) are the radial parts of the J = 3/2 and 
J = 1/2 state wave functions for the electron, respectively. With these definitions, the quantity 
Ck( r i) becomes dimensionless. 

As mentioned above, by the perturbation due to V(r), the states |p 3 / 2 ,/z = ±3/2) are slightly 
mixed with other states ip s and z/> p as: 

\P3/2, P = ±3/2) -± \p3/ 2 ,fJ- = ±3/2) + \ip s ) + | i/j p ) . (8) 

Here 

IVb) = --^jfDiR s (r)\s 1/2 , n') , (9) 

where (i,p!) = (+, 1/2) and (—, —1/2) for fi = 3/2 and —3/2, respectively, and 


IV’p) 


2 \/3* Ry 
5 E p 


Qr 


^2 ( a m)a'P'Rp( r )\Pl/ 2 ,P) , 
±=±1/2 


( 10 ) 
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where (i,/3') = (+, 1),(—, 2 ) for p = 3/2 and —3/2, respectively, and a' = 1,2 for p’ = 1/2 and 
— 1 / 2 , respectively. The definitions of the “±” symbols which appear in the subscripts in Eqs. (J9|) 
and (HOD like D± will be given later in this paper at Eq. (I16D . The factors E s and E p are the 
energy levels of the \si/ 2 ,p') and \pi/ 2 ,p') states {p! = ± 1 / 2 ), measured from the state |_p 3 / 2 , //■) 
(p = ±3/2), respectively, namely, 

E s = E(s 1/2 ,p') - E{p 3/2 , p) , E p = E(p 1/2 ,p') - E(p 3/2 ,p) . (11) 

In Eq. (flOl) . we have neglected the difference between E(p 3 / 2 , 1/2) and E(pi/ 2 ,— 1/2). The func¬ 
tions i? s (r) and R p (r ) are the radial part of the s- and p-wave-state wave functions, respectively. 
The factor a m (m = 1,... ,3) denotes the Pauli matrices, and (a m )ot'i3' ( a '■,$' = 1,2) its ( a',/3') 
component. Note that our overall phase convention of \p\/ 2 ,p') differs from that of Ref. [1J by a 
factor of i. Since 4> s in Eq. © depends on the value of p, we denote the ip s which mixes with 
|p 3 / 2 , p = ±3/2), and \p 3 / 2 ,p = —3/2) as i/ts+ and V’s-, respectively. We also define i^ p+ and ijjp- 
similarly, and use the notation ijj s ± and i/j p ± instead of and ifj p where the distinction is necessary. 
In the subsections just below, we prove Eqs. (JU|) and ffTUl) . 

3.1 Verification of Eq. ([9]) 

In this subsection we verify Eq. Q. To do so, we need the solutions of the Dirac equation for 
an electron in the Coulomb potential, V(r) = — Za/r. Since we are interested only in the non- 
relativistic limit, in this paper we consider only the upper two components of the four-component 
Dirac spinor. 

As a preparation, we first introduce some notations. We define the two-component spinors Xa 
(a = ± 1 / 2 ) as, 


Xl/2 20F (o) ’ X " 1/2 20F (i) ' (12) 

Namely, the angular parts of the states |si/ 2 ,o:) (a = ±1/2) are related to Xa by 

l s l/2) a ) = Xa , (13) 

where we have suppressed the radial part of the wave function. Below we suppress the radial part 
when there may arise no confusion. We also define the vectors e t (i = ±, 0) as 

e+ = -i(l,i,0)/>/2, e_ = i(l,-i,0)/V2 , eb = *(0,0,1). (14) 

They satisfy the relations, 

e+ • e + = e± ■ e_ = 0 , e + • e_ = 1 , e± ■ e 3 = 0 , eb • e*o = — 1 . (15) 

We denote the inner product of a general vector k and e±p as k±p: 

ki = k ■ Si , (* = ±, 0). (16) 
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With these definitions, the quantities ki {i = ±,0) form a spherical tensor of rank 1. Since the 
vectors e) (i = ±, 0 ) are linearly independent of each other, any three dimensional vector k can be 
decomposed in terms of e) as 

k = {k ■ e_)e+ + (k ■ e+)e_ — (k ■ eb)eb 
= k-e+ + k + e- — k^eo . (17) 

Now, by using these notations, the angular part of the \ps/ 2 ,p) (p = ±3/2) states can be 
written as 

\P3/2, p) = VSriiXa , (18) 

where n* = n ■ e) = (r • e))/r with i = ±. The indices i and a in the above equation should be 
understood as (i, a) = (+, 1/2) for p = 3/2, and (i, a) = (—, —1/2) for p = —3/2. 

We are now ready to evaluate the angular integral in the factor (si/ 2 , l/ 2 | 2 D m ro m |p 3 / 2 ,3/2), 
where the integral is to be performed only over the angular variables. The integral reads: 


(si/ 2 , l/2|2T) m n m |p 3 / 2 ,3/2) = 2Z 1 Ci(r 1 ) [ — — (-*) \ sin0e*^ 

' ' J z-y/7r r\ r V o7r 


+ (( r ij-^i) (^ 2 ,^ 2 )), 


(19) 


where dfl = sin 6d6d(f>. The integral is straightforward, and the result is 


( s i/ 2 > l/2|2H m n m |p 3/2 , 3/2) = 2Z 1 C 1 (r l ) { + ( ( rijZl ) (r 2 ,Z 2 ) 


= '‘ + = 7S D+ ■ 


( 20 ) 


Once the above matrix element is determined, the other matrix elements {si/ 2 , p'\2D m n m \pz/ 2 , p) 
for (p, p') = (—3/2, ±1/2), (3/2, —1/2) can also be determined by using the Wigner-Eckart theorem. 
To do so, first by using Eq. (P2D, we write the quantity D m n m as 


DrrJim = D_(e + • n) + D+(e_ • n) - D 0 (e 0 ■ n) 

= D~n + + D + ri- — DqUq . (21) 

Then we write the matrix elements as 

(s 1 / 2l p'\2D m n m \p 3/2 ,p) 

= 2Z?_(si/ 2 , p'\n+ |p 3 / 2 , p) + 2D + (s 1/2 , p\n-\p 3/2l p) - 2D 0 (s 1/2 , p'\n 0 \p 3 / 2 , p) 


= 2i(-l) 3/2 ^ (s 1/2 \\n\\p 3/2 ) 



A/2 1 3/2\ 

\~p' 1 p ) 


+ D+ 



1 3/2\ /1/2 1 3/2\ 1 

-1 P ) ° 0 M / J 

( 22 ) 
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where 


are the Wigner 3j symbols. In the second equality we have used the 


(ji h j's \ 
l mi m2 m3 J 

Wigner-Eckart theorem, 


(j'm'\fk q \jm) = i)jmax—m' [ / k 3 ) (/||/ fc ||j) , (23) 

\—m q ml 

where j max = max(j, j’) and f^ q (q = —k, —k + 1, • • • ,k) are a spherical tensor of rank k. Since we 
already know the value of (si / 2 , l/2|2E> m n m |p 3 / 2 , 3/2) from Eq. (1201) . we can determine the value of 
the double-line matrix element (si/ 2 ||n||p 3 / 2 ) by applying Eq. (f 22l) to the case (//,//) = (3/2,1/2). 
This can be done with the help of the numerical values of the Wigner 3 j symbols as 

2 

(si/ 2 |M|p 3 / 2 ) =i-^= . (24) 

Then, by using the value of (si/ 2 ||n||p 3 / 2 ), the other relevant matrix elements can be determined 
from Eq. (l22ll as 


{ s i/2i l/2|2.D m n m |p 3/ / 2 , 3/2) — , (25) 

( s i/2> — l/2|2D m n m |p 3 /2, +3/2) = (si/ 2 ,+l/2|2ZD m n m |p 3 /2, — 3/2) =0 . (26) 

By combining all the above results, we arrive at 

2 

'y ^ l s i/ 2 > I 1 )( s i/ 2 ; A 4 |2T> m n m |p 3 / 2 , A 1 ) = ~DiXa 1 (27) 

+=± 1/2 ^3 

where i and a should be understood as (i,a) = (+,1/2) for [1 = 3/2, and (i,a) = (—,—1/2) for 
[i = —3/2. From this equation, Eq. (|9|) immediately follows. 

3.2 Verification of Eq. (HOj) 

In this subsection we verify Eq. dm The matrix element we would like to evaluate is, 
(Pi/2,^'|(-3)(3mnn m n n |p 3/ 2,±3/2) = -3ZiC , 2 (ri)(p 1/2 ,// , |(ni • n) 2 |p 3/2 , ±3/2) 

± ((ri,Z{) {r 2 ,Z 2 )\ . (28) 

Note that in our notation, which is actually the notation of Landau-Lifshitz, the states |j>i/ 2 ) m) 
(fi = ±1/2) are related to Xa (a = ±1/2) by 

\Pi/2,v) = -*(<?■ n)Xa , (29) 
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where a = 1/2 and — 1/2 for yu = 1/2 and —1/2, respectively. The overall phase convention of 

|Pi/ 2 ) A 4 ) differs from that of Ref. [1] by a factor of i. 

We first evaluate the factor (pi/ 2 , p'\(Hi • rz) 2 |p 3/2 , ±3/2), where the integral is to be performed 

only over the angular variables. 

To do so, we write the inner product Hi ■ n as 

x y z 
m-n = n lx - + m y - + n lz - 

= nix sin 6 cos </) + n\ y sin 6 sin 4> + niz cos 6 

= i sin 6 f (nix + ini y )e~ l< ^ + (ni x - iniy)e 1 ^ j + ni z cos 9 . (30) 


We now substitute the above expression into the matrix element (pi/ 2 , l/2|(ni • n) 2 |p 3 / 2 ,3/2) 


(Pi/ 2 ,l/ 2 |(ni • n) 2 |p 3 / 2 ,3/2) = j dn^=cos9 (Hi • n ) 2 (-i)^J —sin 9 e l</> 


= j cos sin (nix + ini y )e + (ni x — iniy)e l< ^ J + ni z cos l 



x \l — sin 9e 1 ^ . 
8i r 


(31) 


The angular integral is straightforward, and we are left with: 


(Pi/ 2 . V 2 l(«i ' ^) 2 |^ 3 / 2 ,3/2) 


V 2 _ 

5\/3 


(^la: ± 




(32) 


We now determine the other relevant matrix elements (p\ / 2 , ±l/2|(ni -n) 2 |p 3 / 2 , ±3/2) by using 
the Wigner-Eckart theorem. To do so, we have to rewrite (Hi ■ H) 2 in terms of spherical tensors. 
We construct a rank-2 spherical tensor N 2q (g = —2, —1, • • • , 2) by combining two rank-1 spherical 
tensors n* (i = ±,0). This can be done by using the Clebsch-Gordan coefficients, or equivalently 
the 3 j symbols (see Eq. (107.3) of Ref. [9])- The results are: 


N 2 ,±2 = n± , N 2 ,± i = V2n±n 0 , 


N 2 ,o = -j=n+n- + -y= no 


(33) 


We also use a spherical tensor of rank 0, IVo,o, which we construct according to Eq. (107.4) of 
Ref. [9], 


N 0fi = 2n + n_ - n 2 0 . 


(34) 






We now write {n\ • n) 2 in terms of A/g and A/q: 


(ni • n) 2 = (ni+n_ + ni_n+ — ni 0 n 0 ) 2 

= (ni+) 2 N 2-2 — V2ni + n w N2-\ + 02,0^2,0 ~ V^rai-nioA^,! + N 2 ^ + cq ,0-^0,0 


(35) 


where C 2 ,o and co,o are linear combinations of ni + ni_ and nio^io, whose explicit forms we do not 
need for our purposes here. 

We now apply the Wigner-Eckart theorem. First we expand (n\ ■ ft) 2 in terms of the tensors 

N, 

(Pi/ 2 ,n'\(rii -n) 2 \P3/2,v) = (m+) 2 (pi/ 2^'\N2-2\P3/2, p) - \^n 1 + n w (p 1/2 , p'\N2-i\p 3 / 2l p) 

+ C2,o(Pi/2,p'\ n 2,o\p 3 /2,p) ~ p'\N 2 ,i \p 3 /2, p) 

+ (ni-) 2 {pi/ 2 , p'\N 2 ,2\P3/2, p) + Co ,0 {Pl/ 2 , p' I No,o\p 3 / 2 , p) • (36) 

We use the Wigner-Eckart theorem to rewrite this as 


(pi/ 2 .M'l(«i'«) 2 b 3 / 2 .M)=i 2 (-l) 3/2 -'‘'-|(ni + ) 2 ( 1/2 , 2 3/2 ) -V2n 1+ n w ( 1/2 , 2 3/2 

\-P -2 p J \-p' -1 p , 


I 1/2 2 3/2 \ r- ( 1/2 2 3/2' 

+ c 2 ,o . n - V2ni_m 0 

\-P 0 P \—p 1 P , 


2 , 1/2 2 3/2 


+ (« 1 -) 2 ^ J j(Pl/2\\N 2 \\P3/2) 

+ (-l) s / 2 ->%„( 1/2 ° 3/2 ) (pvjIIJVoIIk/s) ■ 
\-P 0 p J 


(37) 


We can fix the value of (P 1 / 2 II-W 2 IIP 3 / 2 ) from Eq. (1321) as 


(P 1 / 2 IIJV 2 IIP 3 / 2 ) = - 


2y/ld 

5y/3 


(38) 
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By using this value, we can determine the following matrix elements: 


(Pi/2,-l/2|(ni -n) 2 |p 3 /2,3/2) = i 2 (-l) 3/2+1/2 (m+) 2 (Pi/ 2 1|^a|IP 3/2 


2\/2 / __ 2 


5^3 


(«i+)' 


(39) 


(Pl/2,l/2|(m -n) 2 |p 3 /2,-3/2) =i 2 (-l) 3/2 1/2 (m-) 2 f2 ^ 9 j (Pl/2l|W 2 ||p 3 /2} 


2\/2 / 2 


5\/3 


(ni-) 2 , 


(Pi/ 2 ) —1/21(ni • n) 2 |p 3/2 ,-3/2) = i 2 (-l) 3/2+1/2 (-\/2ni_nio) ^ -Z/ij (Pl/ ‘ 


5^3 


ni_nio 


(40) 

211-^21 | P 3 / 2 ) 

(41) 


The results Eqs. (f32l) . (1391) . (1401) . (14TT) can be compactly summarized as 


(Pi/ 2 >^'l(ni ■ n) 2 \P3/2, t*) 


2 i 

5\/3 


nij(a • ni) a ipi , 


(42) 


where (i,(3') = (+, 1),(—,2) for p = 3/2 and —3/2, respectively, and a! = 1,2 for p' = 1/2 and 
— 1/2, respectively. By substituting the above equation into Eq. (f28l) . we obtain 


{Pl/2> P |( 3)Q mrl 7r m 7l n |p3/2, i3/2) — Qim{&m)a'P' j 

where (i,/3') = (+, 1),(—,2) for p = 3/2 and —3/2, respectively, and a' = 1,2 for p! 
— 1/2, respectively. It follows that, 


(43) 
1/2 and 


Y \Pi/2,p'){Pi/2,fJ’'\{-3)Q m nn m n n \p 3 /2 1 ±3/2) = -^-Q im Y (°m)a'p{a ■ n)xn' , (44) 

,u'=±l/2 /i'=±l/2 


where (/,/?') = (+, 1),(—,2) for p = 3/2 and —3/2, respectively, and a' = 1,2 for p' = 1/2 and 
—1/2, respectively. 


4 Evaluation of Anapole Moment 

We are now ready to calculate the anapole moment of the four-atom molecule of Fig. [TJ 
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What we are interested in is the expectation value of a sandwiched by the state represented 
by the right-hand side of Eq. (0). The non-trivial lowest-order contribution comes from the terms, 

(a) = i^ s+ \a\^ p+ ) + (c. c.) (45) 

for n = 3/2, and 


(a) = (V’s-lalV’p-) + (c. c.) , (46) 

for fi = —3/2, where (c. c.) stands for the complex conjugate. We need a proof for the above 
statement (i.e. the statement that the first-order perturbation and the rest of the second-order 
perturbation both vanish), but at this moment, we admit Eqs. (|45l) and (1461) . and calculate the 
right-hand sides of these equations. 

The contribution (a sp m) from the “spin current” to (a) is, 


(a. 


spin/ 


en 

2m, 


■ j d 3 r r 2 x (V4<?V , p)) + ( c - c 0 ■ 


(47) 


We represent the integrand in terms of the spherical coordinates, (r, #,</>): 


(aspin) = -7T— [ r 4 drdn\e r + C ° S ° (^o^p)- ]—-d^laei/jp 

2m e J { \ r r sin 6 r sin u 


( 11 

+ eg -d r (ipl(T^ipp) - (ipta^p) -\ - ^d^l<j r ipp) 

\ r rsmU 


+ ( d r (i>\crgi/jp) + \ipla e ^p) - -do^G^) ) \ + (c. c.) , (48) 


where the vectors e r , eg, are defined as 


^ e r \ 

( — / 

dr/ 

dr 

dr 

\ 


^ sin 9 cos 4> 

sin 9 sin </> 

cos 9 \ 

^ e x ^ 

eg = 

dr / 

del 

dr 

dd 


= 

cos 9 cos (/> 

cos 9 sin </> 

— sin 9 J 

e y 

\ e J 

{£*/ 

dr 

d<j> 

) 


\ — sin (j) 

COS (j) 

o / 

W 


(49) 


where e x ,e y ,e z are the unit vectors in the x,y,z directions in the Cartesian coordinates, respec¬ 
tively. The matrices G r ,Gg, a a are the 2x2 matrices defined in such a way that the identity below 
holds: 


g — & x e x T T g z e z — G r e r T Ggeg T G^e^ , (50) 

where cr x ,y,z are the Pauli matrices. The explicit forms of o r q $ are, 


/ cos 6 e ^ smd 
ye 1 ^ sin 9 — cos 9 


( — sin 9 e ^ cos 9 
cos 9 sin 9 


O(j) 



(51) 
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We can greatly simplify Eq. (148p by integrating-by-parts those terms which have d r , do , or 
d 0 . First, integrating-by-parts those terms with d r is equivalent to replacing the operator d r with 
a factor of (—4/r): 



We now integrate-by-parts the terms with <9g. By noting that doe r = eg and dge^ = 0, we obtain: 

(«s P in) = r 4 drdn^e r + eg (^(VflMp) + ^^WK^p)) 

+ e </> ) + (c- c.) . (53) 

\ r r sm 6 J J 

Finally, by integration-by-parts the terms with d by noting that d^er = sin0e^ and <9<^eg = 
cos 0 e^, we are left with: 


(aspin) = - J r 4 drdf! jeg + e^ j + (c. c.) . (54) 


The integral with respect to r is now easy: 


(aspin) = r{s,p ) J dQ^eod’la^p) angular - e ( / ) (V’J'7g^p)angular| + (c. C.) , (55) 


where the subscripts “angular” means that only the angular part should be considered, forgetting 
about the radial part. The factor r(s,p ) is defined in the unnumbered equation just above Eq. (15) 
of Ref. [I] , 


r 

r(s,p) = 

Jo 


dr r 3 R s (r)R p (r ) 


where R s (r ) and R p (r) are the radial parts of ip s and respectively. 
The contribution (a or b) from the usual “orbital current” to (a) is, 


(56) 


(Sorb) = J d 3 r r 2 j + + (c. c.) . 


(57) 












The integral can be written in the spherical coordinates as: 


(®orb) — 


e7T 

2 m e 


r 4 drdtl jv’K - *) (e r d r ip p ± eg^dgip p ± e 0 r 
+ i (e r d ri’l ± eg^dgipl ± e ^—t/p j + (c. c.) . 


(58) 


We can use d c pip s = dgi/j s = 0 to drop two terms in the second line. Upon integrating-by-parts those 
terms with dg or d^, we find that such terms cancel with each other. We are then left with: 


(®orb) — i 


eir 
2 m e 


r 4 drdO, e r < ipl(d r ip p ) — (d r ipl)ip p > + (c. c.) . 


(59) 


The r-integral is now easy: 


(aorb) = i^-r^s^) I dn e r {^ p ) lar + (c. c.) , 


2m, 


angular 


(60) 


where r\(s,p) is defined as 



Before going further, it is convenient to calculate the inner products (sc 2) M = ±l/2|e r |p 1 / 2) pt' = 
± 1 / 2 ), and (si/ 2 ,// = ±1/2| (ega,/, — e$ 0 g)\pi/ 2 , = ± 1 / 2 ), where the integral is to be performed 

only over the angular variables, for all the possible combinations of pi and pi'. 

( s i/ 2 i d = ±l/2|e r |pi/ 2 , fi' = ±1/2) are evaluated to be, by using the explicit solutions for the 
Dirac equation for a particle in the Coulomb potential, 

( s i/ 2 j +l/ 2 |e r |p 1 / 2 , ± 1 / 2 ) = — * J dD —^1 0 ^ e r 

f 1 

= —i I dQ— cos 0 (sin 6 cos 4>e x ± sin 6 sin (j)e y ± cos 9e z ) 

= ~~x e z ■ ( 62 ) 
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Similarly, 


(si/ 2 ,+l/ 2 |e r |pi/ 2 , — 1 / 2 ) = J dQ—sinde ^ (sin 9 cos (f>e x + sin 0 sin 0 e y + 


cos 9e z ) = 


(si/2,-l/ 2 l e rbi/2>+ 1 / 2 ) = J dn^-smde lll> (;, 


= / dn—sin0e*^(sin0 cos (pe x + sin#sin</>e„ + cos 9e z 
47r 


( s 1 / 2 5 —l/ 2 | e rbl/ 2 )- 1 / 2 ) — y dhl— cos #(sin 0 cos <pe x + sin 0 sin (j)e y + cos 


g( e x * e y) 

(63) 


2 (®a; 


(64) 


+ 3 e , • 


(65) 


Similarly, the other elements (si/ 2 > + = ± 1/2| (e^cr^ — e<j>Go)\P\/ 2 i p' = + 1/ 2 ) can be computed 
to be: 


2 

( s i/ 2 > + 1/ 2 \( e e^<j> ~ e <j>^e)\Pi/ 2 i + 1 / 2 ) = +g e2 > 

( 66 ) 

2 

( s i/ 2 ) +1 / 2 ( e 6&<j) ~ e <f> a e)\Pi/ 2 i —I/ 2 ) = +^(e x - ie y ) , 

(67) 

2 

( s i/ 2 ) —l/ 2 (eecr^ — e ( j > ag)\pi/ 2 , + 1 / 2 ) = + 3 ( 03 ; + * e y ) j 

( 68 ) 

2 

(si/ 2 ,-l/ 2 |(e e <T^-e^<T fl )|p 1 / 2 ,- 1 / 2 ) =--e* . 

(69) 


From Eqs. (162H691) . the sum of the two contributions (a S pin) (Eq. (l55l) f and (a or b) (Eq. (1601) 1 
can be expressed in a compact form: 

(®) = (®spin) 4“ (fl-orb) 

= ( r l( S >P) - 4r ( S >P)) J dQ e r (V4V’p)angular + ( C ‘ C 0 ' ( 70 ) 

We now evaluate the matrix element. First, the contribution from |si/ 2 >l/ 2 ) and bi/ 2 >l/ 2 ) 
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to (?/> s+ |a|V’p+) + (c. c.) is, 


i I C?f^('0s-|- , 0p+)angular®r (c. 


|a,1/2),|p,1/2) 


= ‘T- { (■- *)-<*)) } 

X {(^|^C 2 (r,)n l2 =!<=t±^) + ((*,*) -> (r 2 ,Z 2 ))l + (c. c.) 


= e. 


4 (Ry) s 

: 15 E,Er, 


Z\ Z-2 


n 3 • ni x n 2 


Ci(r 1 )C 2 (r 2 )n 2 z ~ C 2 (ri)Ci(r 2 )ni z 


(71) 


where we have used n 3 = e z . If we further use J • n 3 = 3/2 which holds for the case fi = 3/2, we 
get 


dfi(V> + 


S+ 


V , p+) angular4“ (c> C.) 


|s,i/2),|p,i/2> 


= e 2 (J • n 3 ) 


8 (Ry) s 

45 E.Er, 


Z\Z 2 


n 3 • ni x n 2 


x |ci(ri)C , 2 (r 2 )n 22 - C' 2 (ri)Ci(r 2 )ni 2 | . 

(72) 
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The contribution from |si/ 2 ,l/2) and |_p!/ 2 , —1/2) to ('il> s+ \a\'ip p+ ) + (c. c.) is, 


i I C?f^('05-|-V ; pH-)angular^r (c. C.) 


| s , 1 / 2 ),|p,-i/ 2 ) 


= i-^(e x - iey) | ^ - -^=~-ZiCi(ri) l(jllx ^ + ^(n, Zi) -> (r 2 , Z 2 ) 


x < I — ^-^r-ZiC , 2 (ri)(ni 3; + iniy) ^ nLt + mi ^ ) + ( ( ri , Zi) ->• (r 2 , Z 2 ) ) }> + (c. c.) 


5 E„ 


V2 


4(Ry) 2 
15 E s Ep 


ZiZ 2 


77-3 • \ni x n 2 


+ 


4(Ry) 2 

15-£/ s -£/p 


C'i(ri)C 2 (r 2 )(n 23 ; e a; + n 2y e y ) - C 2 (n)Ci(r 2 )(m x e x + m^) 
Z 2 C'i(n)C' 2 (n)(m xn 3 )|ni Xn 3 | 2 
+ ZiZ 2 Ci(ri)C 2 (r 2 )(n 2 x n 3 ) ^(ni x n 3 ) • (n 2 x n 3 ) 

+ Z 1 Z 2 C'i(r- 2 )C' 2 (ri)(ni x n 3 ) ^(ni x n 3 ) • (n 2 x n 3 ) 

+ ZlCi{r 2 )C 2 (r 2 )(n 2 x n 3 )|n 2 x n 3 | 2 . (73) 

By substituting the sum of Eqs. ([7T]) and (I73|) into Eq. (170]) . we obtain, for fi = 3/2, 


(a) = — fn(s,p ) - 4r'(s,p) N ) ZiZ 2 (J- n 3 ) 

m e V / 45 E S E P 


n 3 • ni x n 2 


+ — (n(s,p) - 4r(s,p) N ) (/• n 3 ) 

m e V / 45E s E p 


Ci(n)C 2 (r 2 )n 2 - C 2 (ri)Ci(r 2 )ni 
^C 1 (r 1 )C 2 (r 1 )(n 1 x n 3 )|n! x n 3 | 2 


+ ZiZ 2 Ci(ri)C 2 (r 2 )(n 2 x n 3 ) ( (ni x n 3 ) • (n 2 x n 3 ) 


+ Z 1 Z 2 C 1 (r 2 )C 2 (r 1 )(n 1 x n 3 ) ( (ni x n 3 ) • (n 2 x n 3 ) 


+ Z|Ci(r 2 )C 2 (r 2 )(n 2 x n 3 )|n 2 x n 3 | 


(74) 


where we have used J ■ n 3 = 3/2 which holds for the case p = 3/2. The terms in the second 
and the subsequent lines are all perpendicular to the axis of the diatomic molecule and 

vanish by averaging the orientation of the molecule with this axis fixed. As long as the rotational 
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symmetry around the 2 -axis is a good symmetry, this averaging is legitimate since there is a 
quantum-mechanical uncertainty relation between J z and the azimuthal angle (j) which describes 
the orientation of the molecule with the A3A4 axis fixed at the 2 -axis. Thus we obtain Eq. ((2]). 
The same conclusion can also be obtained for fi = —3/2. 

In the above discussion, we have neglected the possible contribution from the terms at first 
order of perturbation. Below we show that this can be justified. At first order of perturbation, the 
possible contribution to (a) comes from 

(a) = {ips±\a\p 3 / 2 , I* = ±3/2) + (c. c.) (75) 

(The matrix element of a between \i/j p ±) and |j> 3 / 2 ) vanishes because of parity). To evaluate Eq. (T75L 
we can use 

{ a ) = (dspin) ± ( a orb) (76) 


(®spin) — ^ r(s,J)) J dCl ('i/ , Ju'0'0p) an g U i ar ®</>(V4 , ' r 0V , p)angular^ ± (c. C.) , 


( S orb) = J dn Gr angular + ( C ‘ C 0 ’ 


(77) 

(78) 


where ip p should be understood as \p 3 i 21 P = ±3/2). The integrals which can appear in Eqs. (1771) 
and (17H1) are: 


( s i/ 2 , +l/2|e r |;?3/2, ±3/2) — —j={e x -\-ie y ) , (79) 

( s l/2> ~~ ^-/‘A e r\P 3 / 2 i 3/2) = -J=(e x — i&y) , (80) 

( s i/2) +l/2|e r |p3/2, — 3/2) = (sp/2, — l/2|e r |p 3 /2, ±3/2) = 0 , (81) 

and 

, +l/2|(e0<70 — e^cre)|p3/2, ±3/2) = —^=(e x + ie y ) , (82) 

( s i/2> — l/2|(e£K7^, — e<j>crg)\p 3 / 2 , —3/2) = ~^={e x - ie y ) , (83) 

( s i/ 2 j ±l/2|(eg(T0 — G ( / > CT 0 )\p 3 / 2 , —3/2) = (s 1/2 , —1/2| (ega^ - e^cr 0 )|j> 3/2 , +3/2) =0 . (84) 

From Eqs. (179]) - ([841l . we hnd that the vector ('ip s ±\a\P 3 / 2 i P = ±3/2) + (c. c.) is zero or per¬ 
pendicular to 713, namely, the axis of the diatomic molecule A 3 A 4 . In this case, this vector vanishes 
when averaged over the orientation of the molecule with 713 being fixed. Therefore, to first order 
of perturbation, there is no contribution to the anapole. The same comment also applies to the 
contribution which comes from second order of perturbation like the term 

(«) = (V ; s±I«It’ 3/2,^ = ±3/2) + (c. c.) , (85) 
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where |^) is the s-wave state which appears as a second-order correction to the initial wave function 
|p 3 / 2 ). (Another second-order contribution (ijj' p \a\p 3 / 2 , l 1 = ±3/2) + (c. c.), where ip' is the p -wave 
state which appears as a second-order correction to the initial wave function \p 3 / 2 ), vanishes from 
parity conservation.) In this case, the angular integrals involved are Eqs. (1791) . (l80l) . (1HT1) . (1821) . 
(j83D . and (1841) . which all vanish when averaged over the orientation of the molecule with n 3 being 
fixed. Therefore there is no other second-order contribution to Eq. (J2]) which does not vanish after 
the average over the orientations of the molecule with fixed n 3 as mentioned above. 


5 Discussion 

We have analytically derived the direct relation between the anapole moment and the geometrical 
structure of chiral molecules. A four-atom molecule has been studied as an example since it is the 
simplest molecule having chirality. We have focused on a valence electron which is captured by 
the diatomic molecule A 3 A 4 , and treated effects from the other atoms A\ and A 2 on the electron 
as perturbation. We have computed the corrections to the wave function of the electron in the 
\P 3 / 2 ,P = ±3/2) state by the Coulomb interactions from the atoms A\ and A 2 , and calculated 
the anapole moment. Though this method was first developed in Ref. [I], this method has been 
discussed in much more detail in this paper to be accessible to many scientists over wider regions 
and their final result has been corrected. The electron terms for n-atom molecules (n < 3) reserve 
the symmetry of molecules. This is not the case for n > 4 and needs complicated processes SllSl- 
On the other hand, the procedures developed in this paper have given the direct relation between 
the PV interactions and the geometrical structure without wandering into complicated ab-initio 
MO calculations. Of course, these two approaches are complementary and we need some bridge 
between the two approaches. One point which may need improvement in our approach is that our 
approach might seem to depend on the peculiar initial state |p 3 / 2 ). We may generalize it to the 
cases where the initial state is \p\/ 2 ) or \si/ 2 ) state. Starting from |si/ 2 ) state in place of |p 3 / 2 ), we 
obtain 


(a) = 


2 eir 
15 m e 


ri{s,p ) - 4r(s,p) 


±"(^ 1 / 2 ) - E(p 3 / 2 ,3/2) E(si/ 2 ) — E(p 3 /2, —1/2) J E(si/ 2 ) — E{pi/ 2 ) 


(Ry) 5 


x Z\Z 2 


n 3 • ni x n 2 


C\ (n) C 2 (r 2 )n 2 - C 2 (ri)Ci(r 2 )ni 


4e7r 

15m e 


h{s,p) + 2 r(s,p) 


_ 1 _(Ry) 2 _ 

-^( s i/2i 1/2) — E(p 1 / 2 , —1/2) E(si/ 2 , 1/2) — E(p 3 / 2 , 1/2) 


x Z\Z 2 


n 3 • ni x n 2 


C\{ri)C 2 {r 2 )n 2 - C 2 (ri)Ci(r 2 )ni 


( 86 ) 


where C'i(rj) = o 0 dr r 2 R Sl/ 2 (r)R Pl/ 2 (r)[(n/r 2 ) 6 {r - r;) + ( r/r 2 )0{n - r)], where R Sl/ 2 {r ) and 
R P 1/2 ( r ) are radial wave functions of the initial si / 2 state and the intermediate p \/ 2 state, 
respectively. f(s,p ) and ri(s,p) are defined as f(s,p) = J 0 °° drr 3 R Sl/ 2 (r)R P 3 / 2 (r) and r\(s,p) = 
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/ 0 °° drr 4 (i ? Sl/2 (r)Rp 3/2 (r) — R' Si/ 2 (r)R P 3 / 2 (r)), where R P3/2 {r) is the radial wave function of the p 5 / 2 
state. Instead, if we start from \pi/ 2 ) state, we obtain 


2p7T 

(a) = +- 7 — ( n{s,p) - 4 r(s,p) 
15 m e 


E{pi/ 2 ) — E(p 3 / 2 , 3/2) E(pi/ 2 ) - E(p 3 / 2 , — 1 / 2 ) J E(py 2 ) — E(sy 2 ) 


(Ry) s 


x Z\Z 2 


n 3 • ni x n 2 


Ci{ri)C 2 (r 2 )n 2 - C 2 (ri)Ci(r 2 )ni 


4e7r 

+ 15^1 fl ( s ’^) + 2f ( s >P) 


(Ry) s 


E(pi/ 2l 1/2) — £( Si/2 , —1/2) E(pi/ 2 , 1 / 2 ) — E(p 3 / 2 , — 1 / 2 ) 


x Z\Z 2 


n 3 • ni x ra 2 


Ci(n)C 2 ( 72)^2 - C 2 (ri)C'i(r 2 )ni 


(87) 


This is of a similar form to Eq. (| 86 l) . 

The second point which needs discussions is to apply this approach to realistic molecules. This 
is discussed in Appendix. 

This approach has many applications. One of them is the evaluation of the energy differences 
between optical isomers. The PV potential due to the weak neutral boson exchange is given by 

mmm 

V PV = ^QwWi ■ P^ (3) (r M )} + , (88) 

~r\/ ZTTlp 

c a,i 

where Qyy is the weak charge of the a-th atom in a molecule: 

Qw = ( 1-4 sin 2 6 w )Z a - N a , (89) 

with the atomic number Z a and the neutron number N a of the a-th atom, where a runs over 
the atoms composing of the molecule and i labels the electrons. The precedent analytical methods 
have not discussed the direct relation between the energy differences and geometrical structures. 
The model set-up is the same as in the present paper and we may replace the anapole moment by 
Eq. (1S5D . The detail will be discussed in a separate form |T3]. 
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A Candidate molecules 


The present paper discusses that chiral molecules having an unpaired electron would be good sys¬ 
tems for the measurement of the anapole moment, and the analytical derivation of the expectation 
value of the anapole moment Hamiltonian (Eqs. (47) and (57)) for four-atom molecules, which 
are the smallest molecules that have chirality, is given. In order to further simplify the analytical 
derivation, we assumed that the nature of the unpaired electron is mainly described by atomic 
orbitals attached to A. 4 (Fig.l), and the projection of the angular momentum along the A4 — A3 
axis is a good quantum number. Molecules that satisfy these conditions may be those whose A 4 is 
a heavy atom, in which most of the unpaired electron distributes on A4, and whose A\ and A2 are 
lighter atoms than A3 and A4. 

An example of a four-atom molecule with an unpaired electron is the HNOH radical, in which 
one hydrogen of an amine group is removed from hydroxylamine (NH 2 OH). However, it is not a 
chiral molecule because this radical has a planer equilibrium structure m- We have searched for 
various four-atom molecules similar to the HNOH radical by the use of ab-initio molecular orbital 
calculations, and found several molecules that have a chiral structure. They include the FAsSH 
radical, whose dihedral angle, x, at the equilibrium structure is 61.3 degrees. The potential surface 
of this molecule along the dihedral angle, x> based on MP2/6-31G(d) [15] is shown in Fig.2. The 
potential barrier along x is 354.6 cm” 1 , while the vibrational frequency along x under the harmonic 
approximation is 136.4 cm” 1 . The gross orbital population analysis of a simple ROHF molecular 
orbital calculation indicates that about 82 % of the electron spin locates on the 4 p orbital of the As 
atom. Since the ground state of AsS is 2 n 1/ / 2 ) the anapole moment of this radical may be roughly 
approximated by Eq. ([ST]) . 

Another example is the FPOH radical, which has the dihedral angle, x> of 53.0 degrees. The 
electron spin of this molecule, however, distributes widely over the radical, and therefore it will 
be necessary to include higher order terms that we ignored in this paper to estimate the anapole 
moment of this radical. 

So far a spectroscopic study of the HNOH radical has been reported by matrix isolation 
spectroscopy in solid hydrogen M- No spectroscopic information were reported for other molecules. 
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